Comment on "Ground State Phase Diagram of a Half-Filled One-Dimensional 

Extended Hubbard Model" 
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In Jeckelmann argued that the recently discov- 
ered bond-order-wave (BOW) phase H S H of the 
ID extended Hubbard model does not have a finite 
extent in the {U, V) plane, but exists only on a seg- 
ment of a first-order spin-density-wave-charge-density- 
wave (SDW-CDW) phase boundary. We here present 
quantum Monte Carlo results of higher precision and for 
larger system sizes than previously Using a direct 
finite-size scaling of the BOW correlations, we reconfirm 
that the BOW phase does exist a finite distance away 
from the phase boundary, which hence is a BOW-CDW 
transition curve. We only address the existence of the 
BOW phase and focus on a single value, J7 = 4, of the 
on-site interaction. 

We first determine the critical value Vc of the nearest- 
neighbor interaction. Fig. Q^a) shows the T^-dependence 
of the charge susceptibility Xc{q) at the smallest non- 
zero wave-number, q = 2tt/N, for different system sizes 
N (for the definition of Xc, see 0). The narrowing 
of the peak with increasing N and the convergence of 
the height indicate a charge gap for all V except at 
Vc, i.e., in disagreement with [Ij we find a continuous 
phase transition at [/ = 4. The peak position gives 
Vc = 2.1602 ± 0.0003, which improves considerably on 
the estimate Vc = 2.150 ± 0.010 reported in 0. 

We next choose V = 2.10, where according to there 
should be no long-range BOW. In Fig. ^b) we show the 
corresponding correlation function CB(r) at the longest 
distance, r — N/2, in periodic chains with up to iV = 512 
sites. As a function of 1/A^ for large N, the data scale 
linearly to a value which corresponds to dimerization 
<^ ^ \/C^/2 ~ 0.053 (as defined in Q). It is not clear 
why the density matrix renormalization group (DMRG) 
calculations of (1| failed to detect this rather strong BOW 
order. Dimerization was observed only on the transition 
curve to the CDW phase I*] , where we instead find coex- 
isting critical BOW and CDW fluctuations. The origin 
of this discrepancy at Vc could be that Vc was not deter- 
mined to sufficient accuracy in 0. 

In summary, our improved calculations do not agree 
with the phase diagram presented in jj but reconfirm the 
finite extent of the BOW phase Q as first suggested in 
The advantage of the Monte Carlo method we have used 


FIG. 1: (a) Charge susceptibility at g = 2tt/N for iV = 
32,64,128,256, and 512 (curvature increases with A'^). (b) 
BOW correlation at distance r — N/2 versus inverse system 
size {V — 2.10), along with a linear fit to the large- A'^ data. 


0,13 is that results can be obtained for large periodic sys- 
tems, which are better suited for finite-size scaling than 
the open boundary conditions typically used in DMRG 
calculations on long chains. With the recently improved 
simulation algorithm [sj that we have used here, we hope 
to be able to determine the full extent of the BOW phase 
in the {U, V) plane. 
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